Getting into Shape: Efficient Cooling Scheme for Fermionic Atoms in Optical Lattices 
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We propose an experimental procedure to cool fermionic atoms loaded into an optical lattice. 
The central idea is to spatially divide the system into entropy-rich and -poor regions by shaping 
the confining potential profile. Atoms in regions of high entropy per particle are subsequently 
isolated from the system. We discuss how to experimentally carry out this proposal, and perform a 
quantitative study of its efficiency. We find that the entropy per particle, s, can typically be reduced 
by a factor of 10 such that entropies lower than s/ks ~ 0.2 can be reached. Cooling into highly 
sought-after quantum phases (such as an antiferromagnet) can thus be achieved. We show that this 
procedure is robust against variations of the experimental conditions. 
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Rapid experimental progress in manipulating ultra- 
cold atomic gases has provided physicists with in- 
creased control over quantum many-particle systems 
This was recently evidenced by the observation of a 
Mott-insulating phase of fermionic atoms in a three- 
dimensional optical lattice @, 0]- More complex quan- 
tum phases, such as Neel antifcrromagnets, strongly cor- 
related Fermi liquids or spin liquids in frustrated geome- 
tries, could also be realized using cold atoms. However, 
such phases typically emerge in a temperature regime 
lower than currently achievable. In existing experiments, 
the atomic cloud is pre-cooled by evaporation in a har- 
monic trap and, in a second step, transferred into the pe- 
riodic potential of an optical lattice. Loading the atoms 
into the lattice is ideally performed adiabatically, i.e. 
conserving the entropy of the system. Present experi- 
ments indicate an entropy per particle of s ~ ir 2 T /Tp ss 
1.5 — 2 in the limit of a non-interacting Fermi gas 0, S| 
with T D and Tp the system and Fermi temperatures, and 
ks set to one Q. These values of s are well above the 
onset of interesting correlated phases. Thus, developing 
novel cooling techniques for lattice quantum gases, as we 
propose in this Letter, is a crucial step to demonstrate 
that cold atoms can indeed adequately simulate strongly 
correlated condensed matter systems. 

Cooling atomic gases in optical lattices is the focus of 
an increasing number of studies. For bosons loaded into 
an optical lattice, it was proposed to create entropy-rich 
regions that are later isolated from the rest of the system 
These proposals were inspired by earlier experi- 
ments in which an adiabatic deformation of the external 
trapping potential was used to increase the phase space 
density of Bose gases 0, S] • For fermions in the absence 
of a lattice potential, it was suggested to cool the gas by 
taking advantage of a Feshbach resonance 1,(3. For 
fermions loaded to an optical lattice very few propos- 
als have been put forward. Most of them apply to non- 
interacting Fermi gases [ll| or are based on the use of 



a Bose- Einstein condensate as a heat reservoir [12|, 13 1 . 
However, the possible limitation of entropy reduction due 
to inelastic collisions between bosons and fermions has 
not been addressed yet. 

In this article, we propose an experimentally realis- 
tic procedure to cool two-component fermionic mixtures 
in optical lattices. The key idea is to spatially divide 
the trapped fermionic gas into regions of low and high 
entropy per particle by shaping the trapping potential. 
The two regions are then adiabatically isolated from each 
other and the atoms from the entropy-rich regions are 
disposed of. The remaining atoms have a drastically re- 
duced entropy per particle. In fact, we find that the 
system temperature can be reduced by typically one or- 
der of magnitude while retaining half of the particles. 
In addition, the cooling efficiency remains high over a 
wide range of interatomic coupling strengths, initial par- 
ticle numbers and trap anisotropies. Hence, with this 
method, it should be possible to reach highly anticipated 
quantum phases not yet observed. Such phases include 
the Neel antiferromagnet in a cubic lattice and, perhaps 
even more excitingly, spin liquids or other exotic spin- 
disordered phases in frustrated lattice geometries [14j ] . 
Interestingly, for systems slightly away from half filling, 
we can also reach sufficiently low entropy per particle to 
enter the strongly correlated Fermi liquid regime. Fi- 
nally, our proposal, which relies only on adding a limited 
number of lasers to engineer the trap potential, can be 
well integrated into existing experimental setups. 

Cooling scheme Let us begin with a spin-i mixture 
of fermionic atoms pre-cooled in a dipole trap. As a first 
step, we apply a three-dimensional optical lattice poten- 
tial (Fig. [1] (a)). To allow the atoms to thermalize, the 
loading is done in the presence of a finite but weak in- 
teratomic coupling. We also keep the lattice sufficiently 
shallow for the atoms to redistribute efficiently. As a 
second step, we modulate the entropy distribution by 
creating a potential depression, a dimple, in the middle 
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FIG. 1: (Color online) Cooling scheme. (a) The atoms 
trapped in a parabolic profile are loaded into an optical lat- 
tice, (b) A band insulator is created in a dimple at the center 
of the trap. This core region is isolated from the rest of the 
system, the storage region, by rising potential barriers, (c) 
The storage region is removed from the system, (d) The band 
insulator is relaxed to the desired quantum phase, e.g. a Mott 
insulator by flattening the dimple and turning off or pushing 
outwards the barriers. 
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of the harmonic trap. This dimple must be sufficiently 
deep and narrow for fermionic atoms to accumulate in it 
and form a band insulator [l5j |. The entropy per parti- 
cle in this 'core region' is very small. In contrast, in the 
outer region, called 'storage region', the potential profile 
is kept shallow in order to create a low density liquid 
over a wide volume. Under such conditions and at small 
interaction strength, the entropy per particle in the stor- 
age region is very high. We then separate the core and 
storage regions by slowly rising potential barriers, and 
obtain the potential profile shown on Fig. [1] (b). As a 
third step, we remove the storage region (Fig. Q] (c)). We 
are left with a new effective system characterized by a 
very small entropy per particle (l6| . Finally, as a last 
step, the band insulator is relaxed adiabatically into an 
experimentally relevant phase. For example, if the aim is 
to form a Mott-insulating state, the filling can be lowered 
by slowly flattening the potential in the core region and 
by turning off or pushing outwards the barriers (Fig. [T] 
(d))- 

Shaping the potential profile The above procedure re- 
lies on the ability to add a tailored potential profile on top 
of a lattice potential with amplitude Matticc- To modu- 
late the entropy distribution, the global potential, shown 
in Fig. [2] should realize tight trapping in the core region, 
surrounded by a wide shallow ring in the storage region 
isolated from the core by high potential barriers. To pro- 
duce this profile, we envision to use three elements, (i) a 
shallow harmonic trap (either magnetic or optical), (ii) a 
dimple which confines atoms in a small region around the 
trap symmetry axis and helps to create the band insula- 
tor, and (iii) a cylindrically-symmetric potential barrier 
to isolate entropically poor and rich regions. The dim- 
ple (ii) and potential barrier (iii) are produced by red- 
and blue-detuned laser beams respectively, creating at- 
tractive or repulsive dipole potentials. The dimple has a 
Gaussian profile, while the barrier should rather be a nar- 
row annulus. Experimentally this can be realized cither 
by setting a tightly focused laser beam in rapid rotation, 
or by engineering the beam profile using phase plates or 



FIG. 2: (Color online) Occupation number (dashed line, up- 
per panels), entropy per particle (solid line, upper panels), 
and potential profile 118) (solid line, lower panels) in the pres- 
ence of the dimple and barriers, as a function of the transverse 
(left) and axial (right) coordinate. The potential is offset 
such that V = is at the bottom of the dimple. We chose 
the following experimentally realistic parameters: = 0.5, 
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other diffractive optics 17| . Consequently, in addition to 
the lattice potential, the trapping profile is given by 



armonic H~~ ^dimple am or 

with ^harmonic (r) = V h {x 2 + y 2 + 7 2 z 2 ) / a 2 , 
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where Ufh,d,b} are the potential amplitudes, 7 is a mea- 
sure of the anisotropy of the harmonic trap, u>{d,b} are 
the waists of the gaussian laser beams forming the dim- 
ple and barrier, r& is the radius of the cylindrical barrier, 
and a the lattice spacing. 

Efficiency of the procedure The efficiency of the pro- 
posed cooling scheme can be quantitatively estimated un- 
der the assumption that shaping the potential profile is 
an adiabatic process. Possible deviations from adiabacity 
will be discussed later on. Under the adiabatic assump- 
tion the meaningful quantity is the entropy per particle 
rather than temperature itself. The cooling efficiency de- 
pends on how much entropy per particle is left in the core 
region, sc = Sc/Nc, at the precise moment when the 
increasing barrier height causes the two regions to stop 
exchanging entropy compared to the initial entropy per 
particle, st = St/Nt- The quantities Sq/t an d ^c/t 
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are the entropy and number of atoms in the core (C) 
and total system (T). The described situation is shown 
in Fig. [21 At later times, the core entropy remains un- 
changed as the two regions are now isolated from one 
another preventing the backflow of entropy. 

To determine the efficiency of the cooling scheme, we 
describe the two-component mixture of fermions using a 
Hubbard-type Hamiltonian [19( 



H 
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h.c. 



Here c ia and are the creation and annihilation oper- 
ators of the fermions with a = {f, J,}, J is the hopping 
matrix element, U is the on-site repulsion, fii is the lo- 
cal chemical potential and hur — etc- is the number 
operator on site i. All potential profiles are treated in 
the local density approximation (LDA), i.e. assuming a 
spatially varying chemical potential /i; = fio — V\ . To use 
LDA, local densities and entropies must be obtained for 
the homogeneous system. These quantities are calculated 
using dynamical mean field theory [20| |. In particular, the 
entropy is calculated as in [21 1. 

In the upper panel of Fig. [31 we show, for a three di- 
mensional gas (7 2 = 50), that the final entropy per par- 
ticle in the system core, sc, can be reduced by a factor 
of ten as compared to the initial entropy per particle, st- 
This is done while retaining about half of the particles. 
For s T = 1-95 and N T = 12 • 10 4 atoms, about 5 • 10 4 
atoms are kept. The lower panel of Fig. [3J shows that 
the efficiency of our cooling scheme is very stable against 
variations of the interaction strength and initial particle 
number. Clearly, the procedure is most efficient at small 
values of the interaction strength, but the deviations for 
other interaction strengths are small. Experimentally a 
compromise has to be found between a small value giving 
an optimal gain and the time of thermalization for which 
scattering processes must take place. Finally, we made 
sure that this cooling scheme is efficient both for quasi 
two dimensional (large 7 2 ) and three dimensional (small 
7 2 ) systems. However, in two dimensions, to obtain sim- 
ilar iVc's and maintain the same efficiency, larger radial 
sizes for both the core and storage regions are required 
as less particles can be stacked along the z direction. 

The reduction of the entropy per particle by one or- 
der of magnitude as compared to the current experimen- 
tal situation opens the door to study a wealth of unex- 
plored phenomena in cold atomic systems. As an exam- 
ple, in Ref. 2l|, a (pessimistic) lower bound on the en- 



tropy per particle needed to stabilize antifcrromagnetic 
long-ranged order was estimated to be s ~ 0.2. Using 
our cooling scheme, entropies per particle lower than this 
value can actually be reached in the core region. Starting 
from initial temperatures currently accessible experimen- 
tally, T /Tp ~ 0.15 — 0.2 0, Q, system temperatures of 
T /T F « 0.014 - 0.02 are achieved. 
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FIG. 3: (Color online) Upper panel: entropy per particle 
in the core, sc, and ratio of particles remaining in the core, 
Nc/Nt, after cooling, as a function of the initial total entropy 
per particle, s T , for 7 2 = 50, ±j = 0.5 and Nt = 12 • 10 4 . 
Lower panel: s c as a function of St for different interaction 
strengths and total particle numbers at r y 2 = 50. The pa- 
rameters used to shape the trap are the same as in Fig. [5] 
[23 |. The black arrow in the lower panel indicates the pes- 
simistic estimate of the entropy per particle required to reach 
the antiferromagnet for large interaction strengths [2l| . 



Finally, two remarks are in order. First, higher effi- 
ciencies could be obtained by removing more atoms or 
engineering flatter outer regions that can store more en- 
tropy. Second, the weak dependence of the efficiency on 
the initial entropy (Fig. [3j suggests that this procedure 
can as well be performed several times in a row to reach 
very low temperatures. However, as all changes have to 
be performed slowly, the total time required to cool the 
system will grow with the number of repetitions. 

Removal of storage atoms Having shown how the 
proposed cooling scheme can decrease the entropy in the 
core region, we now adress the fate of the entropy-rich 
part isolated from the core by the potential barrier. If 
storage atoms do not disturb the subsequent experimen- 
tal measurements, they can simply be "pushed outwards" 
by dynamically increasing the barrier radius and raising 
its height [23| to avoid "spilling" the storage atoms into 
the core region. However, in many cases, getting rid of 
these atoms or transfering them to a different hyperfine 
state could be advantageous for later detection. Several 
different removal/transfer schemes may be envisaged de- 
pending on the details of the experimental setup. These 
schemes do not need to be adiabatic as the storage atoms 
are already isolated from the core. 

One possible removal scheme relies on applying a linear 
potential gradient — Fx (which could be due to gravity or 
to an intentionally applied magnetic gradient) and weak- 
ening the shallow trap along the cc-direction. Under the 
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influence of the applied force, storage atoms will undergo 
Bloch oscillations 24} interrupted by Landau-Zener (non- 
adiabatic) transitions to higher bands. These transitions 
can lead to outcoupling of "atom bursts" at multiples of 
the Bloch period Tb = h/Fa (25|. Atoms in the core 
region are confined by the combined dimple/barrier po- 
tentials, and the potential gradient merely shifts the po- 
tential mininum by a small amount. To achieve signifi- 
cant outcoupling rates, one should also significantly lower 
the lattice depth along x. In the weak-binding limit, 
the Landau-Zener formula indeed predicts a transition 
rate r out ~ T~e~ AljZ Hi, where A LZ = maA 2 /4h 2 F, 
and where A is the bandgap which should be as small 
as possible. For instance, for a lattice depth of 0.5 
(A ~ 0.2 E R ), a = 266 nm and F/m ~ 10 m/s 2 , we find 
Tout ~ 10 s _1 for 40 K atoms (Alz ~ 3) and essentially 
zero for 6 Li atoms (Alz ~ 144) Hf|- Another possibility 
to decrease the band gap is to excite the storage atoms to 
a higher Bloch band 24] . In order to leave the core atoms 
untouched, the excitation beams should have a "hollow" 
profile (created using the same techniques as the poten- 
tial barriers) to suppress the transition probability near 
the center of the cloud. 



Conclusion We proposed an efficient scheme to cool 
fermionic atoms confined to optical lattices. This cool- 
ing procedure relies on spatially dividing the trapped 
fermionic gas into regions of low and high entropy per 
particle using a complex potential profile. We find that, 
for a two-component fermionic mixture loaded into a cu- 
bic lattice potential, this scheme reduces the system tem- 
perature by typically one order of magnitude while keep- 
ing approximately half of the atoms. The procedure re- 
mains efficient over a wide range of interatomic coupling 
strengths, initial particle numbers and trap anisotropies. 
This method can be used to cool atoms into highly 
sought-after quantum phases such as the Neel antifer- 
romagnet, spin liquids and strongly correlated Fermi liq- 
uids. 
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Deviation from adiabaticity In an experimental 
setup, one has to find a compromise between changing 
the potential profile slowly (which favors adiabacity) or 
quickly (which subjects the system to external disrup- 
tions only for a short time). To approximate the heating 
induced by non-adiabaticity, we perform time-dependent 
simulations of the cooling and subsequent relaxation into 
a Mott-insulating state within an experimentally realistic 
time. As a measure of the induced heating, we determine 
the energy absorbed by the core region during the pro- 
cess. We use the adaptive time-dependent density ma- 
trix renormalization group method 27, [28} to simulate 
the procedure in a one dimensional fermionic system de- 
scribed by the Hubbard model. The simulations follow 
this sequence: we (i) shape the trap (by increasing the 
dimple amplitude and the barrier height), (ii) relax the 
band insulator (by simultaneously pushing outwards the 
barriers, decreasing the dimple amplitude, and adjusting 
the density by changing the parabolic trapping poten- 
tial), and (iii) tune the interaction strength to its final 
value. We assume a linear variation of each parameter 
with time. For a total procedure time of the order of 500 
4 (700 ms for 40 K atoms in a lattice with Viattice = 8Er), 
the system remains very close to its ground state Hj| and 
the energy absorbed by the system is smaller by more 
than one order of magnitude than the superexchange cou- 
pling 4J 2 /[/. Consequently, the heating induced by the 
non-adiabaticity is small enough not to hinder the effi- 
ciency of our cooling scheme. We expect that for a three 
dimensional system the timescales for the redistribution 
of atoms are even more favorable than for the one dimen- 
sional case simulated here. 
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